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COSMOLOGY IN TERMS OF THE 
DECELERATION PARAMETER. PART I 


Yu.L. Bolotin, V.A. Cherkaskiy, O.A. Lemets 
D.A. Yerokhin and L.G. Zazunov 


’’All of observational cosmology is the search 
for two numbers: Hq and go-” 

Allan Sandage, 1970 


Abstract 


In the early seventies, Alan Sandage defined cosmology as the search for two 
numbers: Hubble parameter Hq and deceleration parameter qq. The first 
of the two basic cosmological parameters (the Hubble parameter) describes 
the linear part of the time dependence of the scale factor. Treating the 
Universe as a dynamical system it is natural to assume that it is non-linear: 
indeed, linearity is nothing more than approximation, while non-linearity 
represents the generic case. It is evident that future models of the Universe 
must take into account different aspects of its evolution. As soon as the 
scale factor is the only dynamical variable, the quantities which determine 
its time dependence must be essentially present in all aspects of the Universe’ 
evolution. Basic characteristics of the cosmological evolution, both static and 
dynamical, can be expressed in terms of the parameters Hq and go- The very 
parameters (and higher time derivatives of the scale factor) enable us to 
construct model-independent kinematics of the cosmological expansion. 

Time dependence of the scale factor reflects main events in history of 
the Universe. Moreover it is the deceleration parameter who dictates the 
expansion rate of the Hubble sphere and determines the dynamics of the ob¬ 
servable galaxy number variation: depending on the sign of the deceleration 
parameter this number either grows (in the case of decelerated expansion), 
or we are going to stay absolutely alone in the cosmos (if the expansion is 
accelerated). 

The intended purpose of the report is reflected in its title — ’’Cosmol¬ 
ogy in terms of the deceleration parameter”. We would like to show that 
practically any aspect of the cosmological evolution is tightly bound to the 
deceleration parameter. 
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In the forthcoming parts we plan to cover the following topics: 

5. Deceleration parameter in different cosmological models 
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• Cosmology with power and hybrid expansion laws 

• Cosmological models with constant deceleration parameter (Berman 
model) 
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6. Kinematic aspects of dynamical forms of the dark energy 

• Kinematics of scalar helds - general description 

• Quintessence 

• K-essence 

• Chaplygin gas 

• Phantom 

• Crossing the phantom divide 

• Polytropic model 

• Quintom cosmology 
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7. Beyond the standard model 

• Modihed gravity 

• f(T) 
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• Models with interaction of the dark components 

• Models with creation of matter 
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Chapter 1 
Introduction 


In the early seventies, Alan Sandage [T] defined Cosmology as the search for 
two numbers: Hq and q^. It seemed too simple and clear: the main term 
in form of the Hubble parameter (HP) determined the expansion rate of 
the Universe and a small correction due to gravity of the matter content was 
responsible for slow down of the expansion. However the situation drastically 
changed at the end of the last century. 

By that time there was a dominating idea that it is the Big Bang model 
complemented with the Inflation theory which is the adequate model of the 
Universe at least at first approximation. However it was in vain hope. The 
cosmological paradigm had to be changed under the force major of the ob¬ 
servations made with ever increasing precision. Since Hubble’s time the 
cosmologists made attempts to measure the deceleration of the expansion 
caused by gravity. Confidence to discover namely the deceleration effect was 
so firm that the corresponding parameter was called the deceleration pa¬ 
rameter (DP). However in 1998 two independent collaborations [21 [3], having 
researched distant supernovae, presented convincing evidence of the fact that 
the expansion of Universe is accelerated. It turned out that the brightness 
decreases in average considerably faster than it was commonly derived from 
the Big Bang model. Such additional dimming means that a given red shift 
corresponds to some additional distance. And it shows that the cosmologi¬ 
cal expansion is accelerated: the Universe expanded in the past slower than 
today. The discovery of the cosmic acceleration was probably one of the 
most importance not only for modern cosmology but also for physics in gen¬ 
eral. The Universe expanding in accelerated way represents the most direct 
demonstration of the fact that our fundamental theories are either incomplete 
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or, even worse, incorrect HE]- 

We would like to stress however that, in spite the great importance of the 
above mentioned discovery, the role of the DP in cosmological dynamics is far 
from being limited by that. The hrst of the two basic cosmological parame¬ 
ters (the Hubble parameter) describes the linear part of the time dependence 
of the scale factor (SF). The linear evolution is dull and monotonous. All 
linear system behave mostly in the same way while every non-linear system is 
’’unhappy in its own manner” (paraphrasing Lev Tolstoi). Linear dynamics 
can be only regular while non-linear one gives much more possibilities. Ex¬ 
istence of bifurcation points in a non-linear system make dynamical regimes 
to change. Presence of local instability brings possibility to exist for chaotic 
regimes. Treating the Universe as a dynamical system it is natural to assume 
that it is non-linear: indeed, linearity is nothing more than approximation, 
while non-linearity represents the generic case. It is evident that future mod¬ 
els of the Universe must take into account different aspects of its evolution. 
However the simplihed model of homogeneous and isotropic Universe based 
on the cosmological principle contains single dynamical variable in form of 
the scale factor and description of its time evolution presents basic task of 
the theory. It is the first terms in the Taylor series for the scale factor a{t) 
in the vicinity of the present time to who are determined by the parameters 
Hq (the linear part) and go (fhe first non-linear correction). 

As soon as the scale factor is the only dynamical variable, the quantities 
which determine its time dependence must be essentially present in all aspects 
of the Universe’ evolution. Basic characteristics of the cosmological evolution, 
both static and dynamical, can be expressed in terms of the parameters Hq 
and go. The very parameters (and higher time derivatives of the scale factor) 
enable us to construct model-independent kinematics of the cosmological 
expansion. 

Time dependence of the scale factor reflects main events in history of 
the Universe. Moreover it is the DP who dictates the expansion rate of the 
Hubble sphere and determines the dynamics of the observable galaxy number 
variation: depending on the sign of the DP this number either grows (in the 
case of decelerated expansion), or we are going to stay absolutely alone in 
the cosmos (if the expansion is accelerated). 

The intended purpose of the report is reflected in its title - ” Cosmology in 
terms of the deceleration parameter”. We would like to show that practically 
any aspect of the cosmological evolution is tightly bound to the DP. We try 
to make the report possibly independent on the current conjecture associated 
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with the discovery of the accelerated expansion of the Universe. The formula 
S = v^t + at^/2 is correct regardless of the DP sign. On the other hand we 
cannot stay aside from the current observational status of the DP. 

In the conclusions of the paper, by commenting about future observational 
values of Hq and go, Sandage wrote: ’’The present discussion is only a prelude 
to the coming decade. If work now in progress is successful, better values 
for both Hq and go (and perhaps even A) should be found, and the 30-year 
dream of choosing between world models on the basis of kinematics alone 
might possibly be realized”. 

A quarter of a century then passed, a new generation of cosmologists 
came, but the world still was in waiting of the ’’correct” model of the Uni¬ 
verse. A hope for happy outcome of the waiting was again associated to 
forthcoming determination of the DP. In 1997 M. Turner wrote [6]: ”At 
the moment, ACDM best accommodates all the observations, but I believe 
the evidence is not yet strong enough to abandon the other CDM models. 
Especially because additional observations will soon be able to decisively 
distinguish between the different models... This observation include DP go. 
LCDM predicts go ~ 1/2, while the other CDM models predict go = 1/2. 
Two groups (The Supernova Cosmology Project and The High z Supernova 
Team) are hoping to determine go to a precision of ±0.2 by using distant 
Type la supernovae {z ~ 0.3 — 0.7) as standard candles. Together, they dis¬ 
covered more than 40 high redshift supernovae last fall and winter and both 
groups should be announcing results soon.” 

Yet two decades almost passed. Both above mentioned groups were 
awarded by the Nobel Prize. However the DP continues to play central 
role both in description of cosmological observations and in theoretical dis¬ 
cussions. And it is in spite of the appearance of the new cosmological favorite 
- the dark energy. As we have mentioned above, just a few decades before 
the cosmology could be defined as a discipline involved in study of only two 
numbers: current Universe’ expansion rate Ho and the DP go. Now we know 
these numbers with reasonable accuracy but it turns out that fate of the Uni¬ 
verse remains still undehned. To predict the fate of the Universe one needs 
not only the current values of the parameters but their time dependence too 
if present. The accelerations are determined by forces and therefore to deter¬ 
mine the Universe’ fate one needs not only kinematics but dynamics of the 
cosmological expansion too. The latter is known to be determined by energy 
content of the Universe, the concept of which had drastically changed in the 
recent time. The exotic dark energy unexpectedly appeared (rather came 
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back) on the scene and now it seemingly governs dynamics of the Universe. 
The dark energy (in form of the cosmological constant or in different dy¬ 
namical forms generated by scalar helds) hnally balanced the energy budget 
and made total density of the energy content in the Universe equal to the 
critical value predicted by the inflation theory. The new energy component 
has negative pressure and as a consequence it causes accelerated expansion of 
the Universe. It affects both the past evolution of the Universe and its future 
history. If the future for the dark energy in form of the cosmological con¬ 
stant is sufficiently clear and it represents monotonic accelerated expansion, 
then for the dynamical forms of the dark energy there is a whole ensemble 
of different possibilities: Big Rip, Big Whimper, Big Decay, Big Crunch, Big 
Brunch, Big Splat etc. 

We stress one more time that the New Cosmology is based on the accel¬ 
erated expansion of the Universe discovered at the end of the last century. 
Although the 2011 Nobel Prize award for the discovery of the accelerated 
expansion of the Universe formally transferred this effect (if this grandious 
phenomenon can be called simply ”an effect”) into the status of physical re¬ 
ality, summarizing the current state of affairs one must honestly admit that 
we now understand very little about the reason why the Universe’ expansion 
is accelerated. The remaining doubts can be overcome by the cross-check 
of the expansion kinematics based on different physical mechanisms. If the 
fundamental conclusion on the accelerated expansion of the Universe turns 
out to be false then it is difficult to expect the consistency of the results ob¬ 
tained in different approaches (expect in the result of very unlikely random 
coincidence). While the most impressive aspect of the cosmological acceler¬ 
ation consists in the fact that the different strategies of its description lead 
to the same conclusion: the Universe entered into the stage of the acceler¬ 
ated expansion. However experience teaches us never relax our vigilance. A. 
Eddington is attributed the saying: do not trust the observations too much 
until they are not conhrmed by theory. This paradoxical notion is not a theo¬ 
retical arrogance, it reflects deeper understanding that not a plain collection 
of facts but their understanding makes the science. General Relativity allows 
the accelerated expansion of the Universe hlled with a negative pressure sub¬ 
stance, but it does not provide yet deeper understanding of this phenomenon. 
Further studies of the kinematics and dynamics of the Universe’ expansion 
will provide long-time work for both the cosmologists-experimentalists (ob¬ 
servers) and theorists. 



Chapter 2 
Background 


2.1 Definitions and simplest relations 

The Hubble parameter and the DP represent simplest cosmographic param¬ 
eters. It is well known [niEH] that the cosmological principle enables us to 
construct metrics of the Universe and make first steps towards interpretation 
of the cosmological observations. Recall that the kinematics is the part of 
mechanics responsible for description of a body’s motion disregarding the 
forces which cause it. In the sense cosmography is just kinematics of cosmo¬ 
logical expansion. In order to build the key characteristic - time dependence 
of the scale factor a{t) - one needs to take equations of motion (the Einstein’s 
equations) and make an assumption about material content of the Universe, 
which allow to construct the energy-momentum tensor. The cosmography 
is efficient because it allows to test any cosmological model which does not 
contradict the cosmological principle. Modihcations of General Relativity 
or introduction of new components (such as dark matter or dark energy) 
evidently change the dependence a{t) but it affects in no way the relations 
between the kinematic characteristics. 

The parameters H{t) and q{t) can naturally be defined making use of 
Taylor series for the scale factor a{t) in vicinity of the present time to: 

a{t) = a (to) + h (to) [t — to] + -a (to) [t — to]^ + • • • (2-1) 

This relation can be rewritten in the form: 

= ^ + Ho[t — to] — y^fo + • • • , (2.2) 

a (^roj ^ 
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Dynamics of Universe is described in frames of General relativity by the 
Einstein equations: 

Rfiu , 

where the energy momentum tensor T^i, describes spatial distribution of 
mass(energy) while components of the curvature tensor and its trace 
R are expressed in terms of metric tensor g^y and its derivatives of hrst 
and second order. The Einstein equations in general are non-linear and 
complicated to solve. The problem is simplihed if one considers mass dis¬ 
tribution with special symmetry properties embedded in the metrics. For 
homogeneous and isotropic Universe described by the Friedman-Lemaitre- 
Robertson-Walker (FLRW) metrics 12.481 the Einstein equations are reduced 
to system of two Friedman equations: 


a?' 


a 

a 


ml, 


(p + 3p), 


(2.3) 

(2.4) 


where p and p are respectively total density and pressure of all compo¬ 
nents present in the Universe at the considered moment of time and Mpi = 
is the reduced Planck mass. This system of equations is incom¬ 
plete: the two equations are not sufficient to completely describe the dynam¬ 
ics of Universe. Lorentz-invariance of the energy-momentum tensor = 0 
leads to the conservation equation 


p + 3H{p + p) = 0. 


(2.5) 


Let us introduce relative density of i-th component of the energy density 

a = Pc = MliH\ D = 

Pc 


Having dehned relative density for curvature 

one obtains the hrst Friedman equation in the following form 




10 






In order to solve the Friedman equations one needs to define the matter 
content in the Universe and construct the state equation for each component. 
In the simplest linear parametrization the state equation takes the form 


Pi = Wipi 


Solution of the Friedman equations for the case of one-component spatially 
fiat {k = 0) Universe with time-independent state equation parameter {w = 
Wi = const) reads 

2 

( f \ 3(l+tn) 

— J , p oc (2.6) 


We normalize the scale factor by the condition a(to) = 1. Such solutions 
exist only in the case w ^ —1. For the Universe dominated by radiation 
(relativistic gas of photons and neutrinos) one has w = 1/3, while tc = 0 for 
the matter-dominated case. Then fl2.6p takes the form 


a{t) oc ( ^ J ! p (X a ^ 

for matter and 

a(t) oc(—J , p (X a ^ 

for radiation. 

The condition p = const (as for cosmological constant) requires w = —1. 
In such a case the Hubble speed remains constant leading to exponential 
growth of the scale factor 

a{t) oc 


Therefore both traditional cosmological components - matter [w = 0) and 
radiation [w = 1/3) - can produce only decelerated expansion of Universe 
(a < 0, g > 0), namely 


q = 


1 for radiation 

1/2 for matter 

—1 for cosmological constant 


(2.7) 


Using the definition of the DP one obtains for the case of spatially fiat Uni¬ 
verse filled with single component with the state equation [p = wp) the 
following result 


q= -{l + 3w). 


( 2 . 8 ) 
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In general, 


fc = 0,±l, p = ^pi, p=^piWi 


and one obtains 



? = Ifl E 


The latter equation can be rewritten in the form 




—1 


1 + 3 


p 

P 


or equivalently 


Q = - 
2 


1 + 3 Wi^li 


k 

^ 2a ?■ 


Let us consider in more detail the model of two-component Universe [T^ . 
Such approximation is sufficient to achieve good accuracy on each stage of 
its evolution. At the present time the two dark components - dark matter 
and dark energy - are considered to be dominating. With this in mind, let 
us represent total density and pressure in the form 


P Pde + Pmi 
P = Pde+ Pm- 


We neglect meanwhile the interaction between the components and as a result 
they separately satisfy the conservation equation fl2.5l) . Let us assume that 
the state equation parameter for each component depends on the scale factor 

Pde = W{a)pde; 

Pm 


With these assumptions 




de — 


1 + 


k 
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These can, formally, be regarded as two equations for fide and flm- Solving 
them we get 


^ _2g-l-3iy fc(l + 3iy) 

™ “ 3(u; - W) ~ 3am^{w-W)' 

^ 2g — 1 — 3ia k{l + 3w) 

™ ^ 3(u; - W) ~ 3a‘^H‘^{w-W)' 

We note that the denominator in these expressions is always positive because 
W < 0 for dark energy. Thus, the fact that flm{de) > 0 imposes an upper 
and lower limit on q, 


i + 3iy 


1 + 


am'^ 


<q< 


1 + 3ry 


1 + 




It is easy to establish relation between the DP and the pressure for the case 
of spatially flat one-component. It follows from the conservation equation 
(I2.5p that 

p w 

p = -. 

3H 1 -I- w 


Using 


one hnds 

For spatially flat one-component Universe with state equation p = wp the 
scale factor is 

2 

a oc 


and therefore 


H 


2 1 
3(1 + w) t 


Using (12.8p one can hnd a simple relation between the current age of the 
Universe and the DP 


^0 


l + q 


(2.9) 
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Let us now give an instructive example of the fact that solution of the 
Friedman equations can be naturally represented in terms of the two param¬ 
eters Hq and go for the case of Universe hlled with non-relativistic matter 
|20j . Let us start with closed Universe [k = +1). The Friedman equations 
can then be represented in the form 



SttG 



a of 
2- + H^ + — = Q. 
a 


( 2 . 10 ) 


Using this equations one can hnd the relations between the current values of 
the Universe’s parameters 


^0 


qo = 


ao(2go - 1)’ 
47iG 

Po- 




( 2 . 11 ) 


Note that in general 


47rG Po + 3po _ Po + 3po _ 3H^ 

* ^3 HS ~ 2p„,^u ■ “ SttG' 

For Universe hlled only with non-relativistic matter one has = 2go- If is 
easy to see that go > 1/2 and flm > 1 as was expected in the closed model. 
Using fl2.1ip one can rewrite the equation for scale factor in the form 




= --l, 

a 


a = 


2go 


i/o(2go-1)3/2- 


( 2 . 12 ) 


It is easy to see that the considered model includes both parameters Ho and 
go. integration of (I2.12p gives 


Substitution 


leads to 




(1 — cost) 


o., . , 

t = — r — sm T 
2^ ^ 


(2.13) 


14 















Because of the relation dt = adr it is evident that the variable r is the 
conformal time. Taking the constants of integration so that a = 0 as t = 0 
(and r = 0) we can see that a = oq at r = tq. Consequently, 


cos To = 


1 - go 

go 


sin To = 


\/2go — 1 
go 


From fl2.13p it follows that the age of Universe in close model is 




arccos 


1 - go \/2go - 1 


go 


go 


The maximum of scale factor reaches at r = vr, 

2go 


— 


/7o(2go-1)3/2' 


As one can see, all parameters of the model can be expressed in terms of 
parameters Hq and go- 

In the case of open Universe the formulae (I2.inp - (I2.12I) transform into 


..2 1 SttG 

2- + H^-\ = 0, 


Hl = 


go = 


1 


ao(l -2go)’ 
A-nG 


(2.14) 


(2.15) 


?>Hl 


Po, 


• 2 _ , 1 « _ 2go 


( 2 . 16 ) 


Now go < 1/2, 0 < VLra < 1. The solution of fl2.16p in the conformal time 
parametrization is 

a = ^(cosh r — 1), t = ^(sinh r — r). 

Current value tq of the conformal time is dehned by the relation 

1 - go 


cosh To = 


go 
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The age of the Universe is 




1 1 — Q'o + — 2go 


In 


go 


go 


In this case again all characteristics of the model are expressed in terms of 
the parameters Hq and go- 


2.2 Classification of models of Universe based 
on the deceleration parameter 

When the rate of expansion never changes, and a is constant, the scaling 
factor is proportional to time t, and the deceleration term is zero. When 
the Hubble term is constant, the deceleration term g is also constant and 
equal to —1, as in the de Sitter and steady-state Universes. In most Uni¬ 
verses the deceleration term changes in time. One can classify models of 
Universe on the basis of time dependence of the two parameters. All models 
can be characterized by whether they expand or contract, and accelerate or 
decelerate: 

(a) H > 0, q > 0: expanding and decelerating 

(b) H > 0, q < 0: expanding and accelerating 

(c) < 0, g > 0: contracting and decelerating 

(d) H < 0, q < 0: contracting and accelerating 

(e) TT > 0, g = 0: expanding, zero deceleration 

(f) H < 0, q = 0: contracting, zero deceleration 

(g) if = 0, g = 0: static 

Of course, generally speaking, both the Hubble parameter and DP can 
change their sign during the evolution. Therefore the evolving Universe can 
transit from one type to another. It is one of the basic tasks of cosmology 
to follow this evolution and clarify its causes. There is little doubt that we 
live in an expanding Universe, and hence only (a), (b), and (e) are possible 


16 



candidates. Evidences in favor of the fact that the expansion is presently ac¬ 
celerating continuously grows in number and therefore the current dynamics 
belongs to type (b). Having hxed the material content we can classify the 
types using the connection between the DP and the spatial geometry. For 
example, in the case of Universe hlled with non-relativistic matter one gets 


1 

g > - ; 
2 

dosed spherical space 

k = +1 

1 

^=2'- 

open flat space 

/c = 0; 

1 

open hyperbolic space 

k = -l 


The classihcation can be also based on the connection between the DP 
and age of the Universe. Foe example, for the Einstein-de Sitter Universe 
where the age is 



we have 

g > ^ ; age< U; 

1 

g = - ; age = t ; 

g < ^ : age> t*. 

Note that sign of the DP determines the difference between the actual 
age of the Universe and the Hubble time. In a decelerating Universe with 
g > 0, the age of the Universe will be less than the Hubble time, because at 
earlier times it was expanding at a faster rate, whereas a Universe that has 
always been accelerating, that is, g < 0 for all time, will have an age that is 
greater than the Hubble time. A Universe that expands at a constant rate, 
g = 0, has an age equal to the Hubble time. 

Finally, if we are interested solely in the expansion regime then in the case 
of constant DP the Universe would exhibit decelerating expansion if g > 0, 
an expansion with constant rate if g = 0, accelerating power-law expansion 
if —1 < g < 0, exponential expansion (also known as de Sitter expansion) if 
g = — 1 and super-exponential expansion if g < — 1. 
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2.3 Deceleration Parameter in the n-dimensional 
Universe 


Following [7], consider an (n + l)-diniensional homogeneons and isotropic 
Lorentzian spacetime with the metric 

ds^ = Qfj^^dx^dx'^ = —dt^ + a?{t)gijdx^dx\ i, j = 1,..., n, (2.17) 

where t is the cosmological (or cosmic) time and Qij is the metric of an 
n-dimensional Riemannian manifold M of constant scalar cnrvatnre charac¬ 
terized by an indicator, k = —1, 0,1, so that M is an n-hyperboloid, the flat 
space M"', or an n-sphere, with the respective metric 

Qijdx'^dx^ = - r— dr"^ + (2.18) 

1 — 

where r > 0 is the radial variable and dVL^_^ denotes the canonical metric 
of the nnit sphere in M”. Inserting the metric fl2.17p - fl2.18p into the 
Einstein eqnations 

(2.19) 

where is the Einstein tensor, G the nniversal gravitational constant, and 
A the cosmological constant, the speed of light is set to nnity, and is the 
energy-momentnm tensor of an ideal cosmological flnid given by 


diag-Ppmj Pmi • • • 1 Pm \) 


( 2 . 20 ) 


with pjn and Pm the t-dependent matter energy density and pressnre, we 
arrive at the Friedmann eqnations 


H 


IOttG 
n(n — 1) 
SttG 
n — 1 


k 

P 2 ’ 


{P + P) + 


k 


( 2 . 21 ) 

( 2 . 22 ) 


in which p,p are the effective energy density and pressnre related to PmiPm 
throngh 

On the other hand, recall that, with fl2.17p and fl2.20p and fl2.23p . the energy- 
conservation law, = 0, takes the form 


Pm + n{pm+Pm)H = 0. 


(2.24) 
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To proceed further, assume that the perfect-fluid cosmological model satishes 
the barotropic equation of state 


Pm ^Pm- 

Inserting (I2.25P into (12.241) . we have 


(2.25) 


Pm + n{l + w)pm- = 0 , 
a 


which can be integrated to yield 


Pm = 


(2.26) 


(2.27) 


where po > 0 is an integration constant [ 8 ]. Using (I2.27p in (I2.23p . we arrive 
at the relation Uli 

p = ( 2 . 28 ) 


From fl2.2ip and fl2.28p . we get the following equation of motion [16] for the 
scale factor a: 


.2 ^ - n ( l +^)+2 ^ y _ p . 

n{n — 1 ) n{n — 1 ) 


(2.29) 


To integrate (I2.29j) . we recall Chebyshev’s theorem [9l|Tn|: For rational 
numbers p, q,r (r 7 ^ 0 ) and nonzero real numbers a, 13, the integral 


x^{a -t- f3x'^y dx 


is elementary if and only if at least one of the quantities 


p + l 
r 


q, 


p + l 
r 




(2.30) 


is an integerli] Consequently, when fc = 0 or A = 0, and w is rational, the 
Chebyshev theorem enables us to know that, for exactly what values of n 
and w, the equation fl2.29p may be integrated. 


^Another way to see the validity of the Chebyshev theorem is to represent the integral 
of concern by a hypergeometric function such that when a quantity in (12.301) is an integer 
the hypergeometric function is reduced into an elementary function. 
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We focus on the spatially flat situation fc = 0, which is known to be most 
relevant for cosmology [m osi ESI [H], and rewrite equation fl2.29p as 

a = ±Vcoa-(i+-)+2 + Aoa2, cq = Aq = (2.31) 

n\n — 1) n[n — 1) 

Then (12.311) reads 

± j a-^ (coa-”(^+"') + Ao) da = t + C. (2.32) 

It is clear that the integral on the left-hand side of (12.321) satishes the integra- 
bility condition stated in the Chebyshev theorem for any n and any rational 

w. 

It turns out that fl2.31l) might be integrated for any real w as well, not 
necessarily rational. To do so we apply a > 0 and get from fl2.3ip the equation 



— In a = ±\/coo ”•(1 +'^) -|- Aq, 
at 

(2.33) 

or equivalently, 




■u = ±\/coe“"'d+'^)“-|-A q, u = \na. 

(2.34) 

Set 

\/coe-’"(i+“’)“ -h Ao = V. 

(2.35) 

Then 

In Co 1 w 2 A ^ 

« = , X , X ln(^^ Ao). 

n(l -|- w) n(l -|- w) 

(2.36) 

Inserting 02.361) into 02.3411. we hnd 



V = T^n(l + w)(n^-Ao), 

(2.37) 


whose integration gives rise to the expressions 

no (l ± ln{l -f w)vot) Aq = 0; 

x/A ;((1 + C'oe^"(^+"')^*)(l - C'oe=F"(i+"')^*)-i, 

Co = ('i^o ~ \/^)('i^o + \/^) Aq > 0; 

\/—Aq tan ^=Fin(l -|- w)\/—Aot + arctan j , Aq < 0 , 

(2.38) 
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where Vq = '^(0). Hence, in terms of n, we obtain the time-dependence of the 
scale factor a: 


n(l+w) 


it) = 


SiiGpo 


n(n 


(2.39) 


We now assume tc > — 1 in the equation of state in our subsequent dis¬ 
cussion. We are interested in solutions satisfying a(0) = 0. 

When A = 0, we combine fl2.38p and fl2.39l) to get 


= AttGpo j (1 + (2.40) 

When A > 0, we similarly obtain 

^ + „)d , (2.41) 

When A < 0 we rewrite (12.391) as 




-A) 


cos 


n(—A) , , /n(n —1) 

-(1 -|- w)t =F arctan \ -—— vq ) . 


2 (n- 1) 


-2A 


If we require a(0) = 0, then (12.421) leads to the conclusion 


(2.42) 


a 




I n(-A) 


■A) 


sm 


2 (n- 1) 


(1 w)t, 


(2.43) 


which gives rise to a periodic Universe so that the scale factor a reaches its 
maximum Om, 


a. 


n(l-\-w) _ 


SiiGpo 

7^’ 


at the times 


(2.44) 




TT 


'2(n- 1) 


{1 + w)y ?7,(—A) 


, fc G Z, 


(2.45) 


and shrinks to zero at the times 


^ io.k 


kn / 2(n — 1) 


{1 + w)y n(—A) 


, k eZ. 


(2.46) 
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Using the equations fl2.40ll2.41ll2.43p . one easily obtains the expressions 
for the DP 



f (1 + “ 1) 


Ao — 0; 



n{l-\-w) 

-1, 

Aq > 0 ; 


q{t) = < 

2 cosh2(,y2(^(i+«'k) 

(2.47) 



-1, 

Aq < 0 . 



[ 2cOs2(,y2(n-l)(l+«'k) 



2.4 Deceleration as a cosmographic parame¬ 
ter 

In the present section we will use an approach to describe the Universe which 
is called a “cosmography” ED- It is entirely based on the cosmological prin¬ 
ciple, stating that the Universe is homogeneous and isotropic on scales larger 
than a hundred megaparsecs. It allows to choose among whole possible va¬ 
riety of models describing the Universe a narrow set of homogeneous and 
isotropic models. The most general space-time metrics agreed with the cos¬ 
mological principle is the FLRW one 

= dt^ — a^{t) I ——i-vT + r‘^{dd‘^ + sin^ 9dip‘^) \ . (2.48) 

( 1 — J 

Here r,6,ip are coordinates of the point which does not participate in any 
motion except the global expansion of the Universe, that is they are the 
comoving coordinates. 

The cosmological principle enables us to build the metrics of the Universe 
and make hrst steps towards interpretation of the cosmological observations. 
Recall that the kinematics is the part of mechanics responsible for description 
of a body’s motion disregarding the forces which cause it. In the sense 
cosmography is just kinematics of cosmological expansion. In order to build 
the key characteristic - time dependence of the scale factor a(t) - one needs to 
take equations of motion (the Einstein’s equations) and make an assumption 
about material content of the Universe, which allow to construct the energy- 
momentum tensor. The cosmography is efficient because it allows to test 
any cosmological model which does not contradict the cosmological principle. 
Modihcations of General Relativity or introduction of new components (such 
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as dark matter or dark energy) evidently change the dependence a(t) but it 
affects in no way the relations between the kinematic characteristics. 

Expansion rate of the Universe is defined by the Hubble parameter H = 
a{t)/a{t) which depends on time. The time dependence is measured by the 
DP q{t). In order to make more detailed description of kinematics of cosmo¬ 
logical expansion it is useful to consider the extended set of the parameters 
which includes higher order time derivatives of the scale factor [221 125112 ^: 


H{t) 

(lit) 

Jit) 

s{t) 

l{t) 


1 da 
a dt ’ 


1 

1 da 

a dt"^ 

a dt 

1 d?a 

1 da 

— 

a dt^ 

a dt 


1 d^a 

1 da 

— 

a dt* 

a dt 


1 d^a 

1 da 

— 

a dt^ 

a dt 



Note that the four latter parameters are dimensionless. In terms of the 
conformal time 



where prime denotes derivative w.r.t. the conformal time. 

We will concentrate our attention on the DP q{t). However one should 
note unique importance of the jerk parameter Jit), in particular for testing 
of the cosmological models. ”It is a striking and slightly puzzling fact that 
almost all current cosmological observations can be summarized by the simple 
statement; The jerk of the Universe equals one” |25] . 

Let us make Taylor expansion of the scale factor in time using the above 
introduced parameters 


a(t) —uq 


1 -|- iLo(^ ~ ^o) 




\qoHoit - to)^ + ^JoH^it - to)^ 


(2.49) 
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Using the same parameters the Taylor expansion of the redshift takes on the 
form 


1 + ;^ 


1 + Ho{t - to) - ^qoHo{t - to)^ + ^joH^it - to)^ 

+ ^’So-hrg(t — to)^ + -^loHo{t — to)^ + O ((t — to)®) 


^ —Ho{t — to) + H^it — to)^ + ■ ■ ■ . 


Derivatives of the lower order parameters can be expressed through the 
higher ones, for instance 


dq 


= j - q{2q + l). 


(iln(l + 2 ; 

Let us give a few useful relations for the DP 

=4 -1; 

q{z) = 


dt 

1 dlnH^ 

2 (iln(l + z) 


, . d\nH , 

=— 

l(a) = - fl + 

dln(all) 


dH 

dt 


= - 1 + 


dH 

da 

H 


q{a) 


- 1 : 


din a ’ 

, dE ^ H 

q — — l + {l + z)E —, E =—] 

dz Hq 


E I{z) = 


1 + q{x) 
1 + a; 


dx. 


For a single-component fluid with density p 


q[a) = -1 - 


da 

2p' 


( 2 . 50 ) 
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The derivatives dH/dz, d^H/dz^, d?H/dz^ and d^H/dz^ can be expressed 
through the DP q and other cosmological parameters 


to 1 0 tq 

1 + 2 
j - 

dz‘^ 

(1 + 2)2 

d^H 

H 

dz^ 

(1 + 2)3 

d^H 

H 

dz^ 

(1 + 2)^ 


(2.51) 


H- 


- Aqj - 3j - s) ; 

(-12g2 - 2V - + 32qj + 23q^j + 7qs + 12j - Af + 8s + l) . 


terms of the scale factor 


1 d^H 


z=0 


^ 2 dz^ 


1 d^H 


z=0 


^ 6 dz^ 


z-' + --- = 


z=0 


= Ho{l + (1 + qo)z + ^{jo - ql)z^ + ^(3^^ + 3gj^ - 4goJo - 3jo - s)^^)- 

2 D 

There is a decomposition for the inverse Hubble parameter 


A 

dz 

dz'^ 



1 dH 1 + g 1 
lP~dz ^ ~l + zH' 


=2 


1 + q 
1 + z 


1 

H 


2 + 4g + 3g^ — j 


(i + zY 


1 


1 

'^0 


1 — (1 + qo)z + 


2 + 4go + 3gQ — jp 
6 


2 ;^ + ... 


Let us cite some useful relations which allow to transform from higher 
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order time derivatives to derivatives w.r.t. red shift; 




= {^ + z)H 


H+{l + z) 


dH 

dz 


d , n d? 

— + + — , 
dz dz^ 

fdH\ 


— = -(l + z)HlH^ + {l + zY[— +{l + z)H 


dt^ 


dY_ 

dt"^ 


\ dz J 


dH ^d^H 


3{1+ zYH^ \H+{1 + z)^ 


= Y + z)H 


,dH 


d? 3 3 d^ 

— -{l + z)H—, 
dH 


H^ + 11(1 + z)H^— + 11(1 + z)H— + il + z) 


dz 


dz 


dH 

dz 


d^H 


\9TTd^H .ndH d^H 

+7(1 + zfH^ + 4(1 + zf — ^ + (1 + zfH^^ 
dz^ dz dz^ dz'^ 


+ {l + zfH^ 


dH 


fdH\‘ 


7H^ + 22H— + 7(1 + zYH — +477 


d_ 

dz 

d^H 


dz 

+ 6{1+ zYH^ \H+{1 + z)^ 


\dz J 

d^ 


dz"^ 


dY_ 

dz"^ 


Derivatives of the Hubble parameter squared w.r.t. the redshift d'^H'^/dz^, 
f = 1, 2, 3,4 in terms of the cosmographic parameters take on the form 


7772 

- (1 

dz 


d^H^ 

2772 

dz'^ 

■(l + ;^)2 

d^H^ 

2772 

dz^ 

■(l + ;^)3 

d*H^ 

2772 

dz'^ 

■(1 + ^)^ 


(l + 2g +j), 

{-qj - s), 

(4gj + 3qs + 3q‘^j - f + As + 1). 


The current values of deceleration and jerk parameters in terms of iV = 
— ln(l + z) read 


go = — 


1 rid(772) 


772 2 dN 


+ 77^ 


v=o 


, _ f 1 d^ ( 772 ) 3 d ( 772 ) 

“ \2772 diV2 ^ 2772 dN ^ 


N =0 
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Time derivatives of the Hubble parameter can be also expressed in terms of 
the cosmographic parameters 

H = -H\l + q), (2.52) 

H=H^{j + 3q + 2), 

H =H‘^[s - Aj - 3q{q + 4) - 6], 

H =H^[l - 5s + 10(g + 2)j + 30(g + 2)q + 24], 

It is easy to see from fl2.52p that accelerating growth of the expansion 
rate H > 0 corresponds to q < —1. 

At last, from the power series expansion of the scalar factor fl2.50p one 
can also express the DP as a power series in time. This time variable can 
be written as a power series in redshift ^ or j/-redshift z/{l + z), yielding 
respectively |2B] 

=?o + {-qo - 2go + Jo) z + ^ (2go + Sq^ + Sq^ - 7qojo - Ajo - Sq) 2 :^ + 0{z^), 

q{y) =qo + (-go - 2go + jo) y + ^ (4go + 8g^ - 7qojo - 2jo - so) y‘^ + 0{y^). 

As can be seen from the relations (I2.5np . the Hubble parameter is connected 
to the DP by the integral relation 

H = Hq exp 

Then it immediately follows that one needs the information on the dynamics 
of cosmological expansion coded in the quantity q{z) in order to reconstruct 
the basic characteristic of the expanding Universe which is H{z). 

2.4.1 Cosmological scalars and the Friedmann equa¬ 
tion 

Dunajski and Gibbons |25] proposed an original way to test the General 
Relativity (GR) and the cosmological models based on it. The procedure 
implies expressing the Friedmann equation in terms of directly measurable 
cosmological scalars constructed out of higher derivatives of the scale factor, 
i.e. cosmographic parameters H,q,j,s,l. In other words, the key idea is 


J [q{z') + l]dln(l + z') 

0 
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to treat the Friedmann equations as one algebraic constraint between the 
scalars. This links the measurement of the cosmological parameters to a test 
of GR, or any of its modihcations (which would lead to different constraints). 
For an example consider a Universe containing the cosmological constant A 
and non-relativistic matter (dust). The Einstein equations reduce to the 
Friedmann equation 

+ fc = ]-pa^ + (2.53) 

o o 

and the conservation equation, which can be used to hnd the dependence of 
density on the scale factor {pa? = M, M = const). Let us now consider a 
system of three equations consisting of fl2.53p and its hrst two time deriva¬ 
tives. We regard this as a system of algebraic equations for the constants 
{k, A, M) which can therefore be expressed as functions of a, a, d, 'a. Take the 
third derivative of fl2.53p and substitute the expressions for {k,A,M). The 
resulting equation does not contain any parameters and can be expressed in 
terms of the cosmological scalars as 

s + 2{q + j) + qj = 0. (2.54) 


This fourth order ODE is equivalent to the Friedmann equation and has an 
advantage that it appears as a constraint on directly measurable quantities. 
Thus it provides the test of the model basing on which the latter equation 
was obtained. 

If only two constants (A, M) are eliminated between fl2.53p and its hrst 
derivative then the second derivative of fl2.53p yields [2^ 

k = a^H\j - 1). (2.55) 


where k is regarded as a parameter. In particular if /c = 0 (SCM) this relation 
reduces to a third order ODE 

3 = 1- (2.56) 

This constant jerk condition is consistent with recent redshift analysis [2^] . 

Let us perform an analogous procedure for the two-component Universe 
hlled with non-relativistic matter with density Mm/o? and radiation with 
density Mr/a^ which do not interact with each other. We represent the hrst 
Friedmann equation in the form 


k Mr, 






+ 


Mr 


SttG 


= 1 . 


(2.57) 
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We then differentiate the latter expression twice w.r.t. time to find 


1 Mm Mr 

2 a2 a3 ’ 




’ I O 

"Cl -|- O" 


Mr 


-a. 


Using the definitions of the cosmographic parameters, one obtains 


+ , = A^3B, (2,58) 

and then 

A = -2j + 6g, B =j - 2q. (2.59) 

In terms of the variables A and B the Friedmann equation fl2.57p takes on 
the form 

^ = {A + B-1)H\ 

or in terms of the cosmographic parameters 


k = a^H‘^{4:q — j — 1). 


Let us check the latter expression for the curvature k in two limiting cases: 
a) a fiat Universe filled only with non-relativistic matter; 


b) a fiat Universe filled only with radiation. 

In the first case (a) one has 5 = 0, g = 1/2, therefore A = 1 and /c = 0. 
Note that j = 1 in this case. 

In the second case (b) one has ^4 = 0, g = 1, therefore B = 1 and fc = 0. 
In this case j = 3. 

The same result can be obtained using the relation fl2.58p and it does 
not require to know value of the deceleration parameter. From the relation 
fl2.58p it follows that in the one-component fiat case A = 1 (a) or i? = 1 (b) 
respectively and therefore j = 1 in the first (a) case and j = 3 in the second 
(b) case. 

Let us now obtain analogue of the equation fl2.54p for the Universe filled 
with non-relativistic matter (dust) and radiation. In this case 


s = -371 




35(1+ 4g). 
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or using the equation 02.591) . one finally gets 


■s - +jQ = 0 . 

It is easy to see that the latter formula reproduces results of the limiting 
cases. In absence of the radiation q = 1/2, j = 1, and therefore s = —7/2. 
This result can be easily verified by direct calculation of the parameter s, 
which can be easily done because a oc in this case. In absence of the 
matter one has instead a oc so q = 1, j = 3 and s = —15. 

2.5 Averaging deceleration parameter 

Since the DP g is a slowly varying quantity (e.g. q = 1/2 for matter- 
dominated case and g = — 1 in the Universe dominated by dark energy in 
form of cosmological constant), then the useful information is contained in 
its time average value, which is very interesting to obtain without integration 
of the equations of motions for the scale factor. Let us see how it is possible 
[29j. For that purpose let us define average value g of this parameter on time 
interval [0 ,fo] with the expression 




1 

to 


to 


0 


q(t)dt. 


Making use of the definition of the DP 


Qit) 


da d / 1\ 

dt \H) ’ 


it is easy to obtain 


Q{to) — —1 + 


toHo 


or 


to — 


Ho\ 

1 + g’ 


(2.60) 


As expected, current age of the Universe is proportional to but the 

proportionality coefficient is solely determined by the average value of the 
DP. It is worth noting that this purely kinematic result depends on curvature 
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of the Universe, nor on number of components filling it, nor on the type of 
gravity theory used. 

The result fl2.60p can be presented in the form 


T 


1 + g’ 


(2.61) 


where T is age of the Universe. As q is of order of unity, it immediately 
follows from (I2.6ip that the Hubble time is s characteristic time scale 
on any stage of Universe’s evolution. Of course the result (I2.6ip coincides 
with that of fl2.9p for the case of constant DP. 

Let us dwell on the result fl2.60p . For the FLRW metric the age of the 
Universe in terms of the redshift ^ is given by the expression 


to — iL, 


-1 


dz 


(1 + 


Taking SCM values for the parameters (Qmo, ^ao), h is easy to show that 


OO 

f _cb__ ^ 

J (1 + + ^AO 

thereby showing that the age parameter for the ACDM model is ifoO ~ 1- 
According to fl2.60p . it is equivalent to g ~ 0. This result can be also obtained 
for one-component liquid with w = —1/3 (the iP-matter [29]), if one recalls 
the result fl2.8p for one-component spatially flat Universe. In such a model 
a{t) = aoHot and d = const, therefore it is called ’’the K-matter coasting 
model”. Another example of coasting model is provided by the Milne empty 
Universe. However, the coasting empty relativistic solution can be achieved 
only in the hyperbolic geometry of the Universe, which clearly contradicts 
the present observations. 

In Fig Em [29| we compare the values of the age parameter Hoto for a large 
set of cosmologies, including ACDM, Einstein - de Sitter, and iP-matter mod¬ 
els. For the SCM (ACDM), the age of the Universe nowadays is exactly the 
same one predicted by the coasting iP-matter model. Naturally, such a fact 
may be just an unexpected coincidence. However, the History of Science has 
already shown that coincidences, mainly in the held of cosmology, deserves a 
special attention. In terms of the average DP the result to = means that 
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Figure 2.1: Age of the Universe in ACDM and other cosmologies. The present 
total age of the Universe for the cosmic concordance model (ACDM) is ex¬ 
actly the same as predicted by the flat coasting iF-matter model, to = 

m- 
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Radiation 



Figure 2.2: The DP as a function of the redshift (the values of z are not in 
scale).The Universe emerges from an inflationary stage at the inhnite past 
(z = 0) and evolves to the inhnite future (z = —1). Abrupt transitions has 
been assumed for any two subsequent stages. In the actual Universe they 
are smooth, but, qualitatively, the result is the same, namely, the DP of the 
Universe seems to be oscillating around the iF-matter solution (g = 0). The 
fate of the Universe is heavily dependent of the next transition in the future. 
The Universe may decelerate (not necessarily mimicking the dust behavior 
suggested in the hgure) [29]. 

the average DP q must be identically zero when it is averaged for a long time 
interval. It should be stressed that it implies equality to zero for the average 
DP which can be achieved through a cascade of accelerating/deceleration 
regimes. An example of such a multiple regime cascade is shown on Fig J2.2[ 
Let us give some remarks on Fig |2.2l which presents a hypothetic history of 
the Universe as a function of the redshift. It starts with the inflation period 
generated by the cosmological constant (q = —1). In order to have the cor¬ 
rect age of Universe (Hq^o = Q = 0), the average value of the DP after 
transition to the accelerated expansion regime until the present day must be 
— 1/2. Concerning the future (at 2 ; < 0), the dilemma is whether the Uni- 
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verse will transit to an accelerating de-Sitter regime {q = —1) as reqnired 
by the LCDM model, or to a decelerating stage which is predicted by some 
scalar held and brane-world scenarios (this qnestion will be discnssed in more 
details in the section ’’Transitive acceleration”). Thns, if g = 0 remains trne 
in the near fntnre, the Universe mnst evolve to a decelerating regime (g > 0). 
Natnrally, this does not mean that snch a transition shonld be the last one. 
We cannot exclnde possibility of a seqnence of transitions with sign change 
of the DP. 

2.6 Energy conditions in terms of the decel¬ 
eration parameter 

Dynamic model-independent constraints on the kinematics of the Universe 
can fnrther be obtained from the so-called energy conditions [201 EB [321 
1331 133] . These conditions, based on qnite general physical principles, im¬ 
pose restrictions on the components of the energy-momentnm tensor In 

choosing a model for the medinm (a model, bnt not the eqnation of state!), 
these conditions can be transformed into ineqnalities restricting the possible 
valnes of pressnre and density of the medinm. In the Friedmann model, the 
medinm is an ideal liqnid, for which 


— (p T P9^u 


where is the 4-velocity of the ideal liqnid of energy density p and pressnre 
p,which can be expressed via the scale factor and its derivatives. 



(2.62) 


In terms of density and pressnre the energy conditions take on the form 


NEC 
WEC => 
SEC ^ 
DEC ^ 


P + P >0, 
P >0, 

p + ?>p >0, 

P >0, 


P + P >0, 
P + P >0, 
—p < p +P- 


Here, NEC, WEC, SEC, and DEC correspond to the zero, weak, strong, 
and dominant energy conditions. Becanse these conditions do not reqnire 
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any definite eqnation of state for the snbstance filling the Universe, they 
impose very simple and model-independent constraints on the behavior of 
the energy density and pressnre. Hence, the energy conditions provide one of 
the possibilities for explaining the evolntion of the Universe on the basis of 
qnite general principles. With expression fl2.62p . the energy conditions can 
be expressed in terms of the scale factor and its derivatives: 


NEC 
WEC => 
SEC ^ 
DEC => 


^- 2 

a 


a 

a 


a 

a 


+ 


> 0 , 

> 0 , 

< 0 , 

>0. 


(2.63) 


In the case of a flat Universe, conditions fl2.63p can be transformed into 
restrictions on the DP q\ 


NEC ^ g > 1, (2.64) 

SEC ^ g > 0, 

DEC g < 2. 

There is no WEC among above conditions becanse it is always satished for 
arbitrary real a{t). 

The conditions fl2.64p considered separately in principle allow a possibil¬ 
ity for both decelerated (g > 0) and accelerated (g < 0) expansion of the 
Universe. The constraints by NEC in (I2.64p have clear sense: as follows from 
the second Friedmann eqnation, the ineqnality p -|- 3p < 0 gives necessary 
condition for accelerated expansion of the Universe, i.e. the accelerated ex¬ 
pansion of the Universe is possible only in presence of components with high 
negative pressnre p < —p/3. The SEC exclndes existence of snch compo¬ 
nents. As a result, g > 0 in this case. At the same time, NEC and DEC 
are compatible with the condition p < —p/3 and therefore they allow the 
regimes with g < 0. 

It worth noting that even before the discovery of the accelerated expansion 
of the Universe in 1997, Visser [35l [36] already concluded, basing on analysis 
of the energy conditions, that current observations suggest that SEC was 


35 












violated sometime between the epoch of galaxy formation and the present. 
This implies that no possible combination of ’’normal” matter is capable of 
htting the observational data. 

Let us make an important remark [37]. In order to connect the energy 
conditions with observations, one often needs hrst to integrate them, and then 
hnd the corresponding constraints on some observational variables, such as 
the distance modulus. Those integral forms can be misleading, and great 
caution is needed when one interprets them physically. 
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Chapter 3 

Distance-Deceleration 
Parameter Relations 

3.1 Different cosmological distances in terms 
of the deceleration parameter 

In cosmology there are many different and equally natural definitions of the 
notion of ’’distance” between two objects or events, whether directly observ¬ 
able or not. Let us consider [38] several examples of distances between two 
objects (events) in cosmology. 

1. The ’’luminosity distance” di. The luminosity distance di of an 
object at redshift z is di = (L/27rF)^/^, where L is the bolometric luminosity 
for a given object and F is the bolometric energy flux received from that 
object. The expression for the luminosity distance in a FLRW Universe is 


/' r z 



/? sinlf ^ r dz 
rrsinn J , 

open 


diiz) = (1 -F z) < 

Z 

ij—i r dz' 

^0 J E(z')^ 

0 

flat 

(3.1) 


HoR 1 E(z') ’ 

closed 



(Here and below we set c = 1). Here R is the (comoving) radius of curva¬ 
ture of the open or closed Universe, E = H/Hq. The relation fl3.ip can be 
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rewritten in terms of the DP. For the spatially flat case 


dz' 


dL{z) = ( 1 + 2 :) J = {l+z)Ho J duexp 


[1 + g(r;)](i[ln(l + n)] 


It is useful to give an expression for the luminosity distance up to terms 
of order of z"^ 


dr = 


Hn 


1 + 


1 - go 


z + 0{z^) 


(3,2) 


where in the spatially flat case 


go — - + SrCj). 


It follows from 03.21) that for small 2 ; the luminosity distance is linearly propor¬ 
tional to the redshift, while the proportionality coefficient equals to inverse 
value of the Hubble constant. For more distant cosmological objects the lu¬ 
minosity distance in the next order depends on current value of the DP go) 
or equivalently on number and type of the components Filing the Universe. 
Expression for the luminosity distance in the next order in the red shift can 
be presented in the form 


di^z) 


cz 



1 + ^(1 - go)^ - ^(1 - go - 3go + 

2 0 

(2 — 2go — 15go — ISgg + 5jo + lOjogo + ■So)2:^ + 0{z^) 


(3.3) 


As was expected, the latter decomposition contains the cosmographic param¬ 
eters dehned in terms of the higher order time derivatives of the scale factor 
{j (X d?a/dt^). 

As an example we calculate the luminosity distance in the Universe Filed 
with non-relativistic matter. Let us start with the Einstein-de Sitter model 
(A; = 0). We represent the spatially Fat case in the following form fl3.ip 


to 

dL = aori{l + z), ri = 



(3.4) 
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In the considered case 


ri = 


3to 

ao 


1 - 


1 / 3 ' 


1 I 

1 + z = — = 

ai 


^ I 

tj 


\ 2/3 


thus 


ri = 


3to 

ao 


(1 + ^) 


- 1/2 


aoHo 


(1 + ^) 


- 1/2 


We took into account here that in the Einstein-de Sitter model hgto 
For the luminosity distance one ultimately obtains 

dh = aon(l + z) = ^[l + z- \/l + z'^ . 


= 2/3. 


2. The’’photon flux distance” dp. Introduction of this kind of distance 
is dictated by the fact, that it is often technologically easier to count the 
photon flux (photons/sec) than it is to bolometrically measure total energy 
flux (power) deposited in the detector. Such transition leads to additional 
factor (1 + z)~^^‘^, thus 

= (i + lyn- 

3. The ’’photon count distance” dp is related to the total number of 
photons absorbed without regard to the rate at which they arrive. Thus the 
photon count distance contains one extra factor of (1 + z)~^ as compared to 
the (power-based) luminosity distance 


dp{z) = 


dr 


1 + ;^ 


4. The ’’angular diameter distance” dA at last is known to be related 
to the luminosity distance in the following way 


dA{z) 


dp 

{1 + zy 


Let us cite expressions for the above considered types of cosmological distance 
up to terms of order z^ [39] : 
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The luminosity distance 


dr — 




+ 2 ;^ • 

+ • 

7joq? 


z + z^ 

1 + ^ 

12 ^ 24 


% 
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^ 5jogo _ ^ 
12 ^ 12 8 


20 40 12 


^0 I 

120 ^ 20 


M , fo\ 

8 24/ 

lljogo ^ ‘27ql 


12 


40 


0 ^ iM + M 


11s, 


0 


120 


+ z" •-h — 

130 30 


1 , 7jo 19j2 , 19/, 


72 


+ 


goSo \ 
8 / 
rr^ 

720 ^ 720 


+ 


% 13jogo 
30 ^ 9 


7jogo 7/ogo _ ^ 133jogg _ 13gg 
18 ^ 240 10 ^ 48 6 

7joqo _ 133gg _ 2^ 1^ _ YjoSp lOgoSp 7q^\ 

4 48 16 90 144 ^ 48 24 / 


(3.5) 


It should be noted that series expansion of fl3.4p in terms of z reproduces fl3.5l) 
for hxed values of the cosmographical parameters, in particular q = qg = 1/2. 


The photon flux distance 


dr? — 


Hp 

4 




+ 


+ 


2 do 3 f 1 Jo 

" ■ " ' Y + " ' ("m “ ? + 12 + 2 ; 

]_ ,]jo _ }Jdo 7>jpqp _ 7^ 

24 ^ 24 48 12 8 
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5^0 , So 


+ 24,' + 


71 


1920 


47jo ^ _ Jo_ 149go _ Ojogp 4)^ 
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_ 7]pql 2^ I ZM _ ^ + 

8 16 8 80 8 / 

+ ^6 , 45 7 jo _ 1^ I _ 1069go 593jpqp 
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7joqi _ 77qp _ 2lql 593so _ 7ipSp llgoSp 7ggso y 
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(3.6) 
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The photon count distance 


dp — 




+ 2 :^ • 

+ ■ 

_ Jhdl 
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, 2/1 5'o\ I 3 

" + " ■ 1 - 2 - 7 ) + " 

-1 + ^ ^ + ^-^0^0 
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_ , _ 4jogo ^ 
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gogo\ 
8 ) 
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8 15 

1 ^ _ 2 ^ 5/0 mo 
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175jogo 
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25goSo 
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(3.7) 


The angular diameter distance 


dA 


1 


;2 + 





4 / 25 13jo 23go 5jogo 13gg 5gg , 

^ V 12 ^ 24 12 12 8 8 ^ 247^ 

5 /137 137jo jj k 163go 7jogo 137gg 

^ V 60 120 ^ 12 120 ^ 60 4 40 

_ 7jogo I _ Zf^ _ I 

8 8 8 40 8 7 


+ / / 49 79jo 31jg 31/0 mp 71go 163jogo 

^ V 20 ^ 40 72 ^ 720 ^ 720 20 36 

7jggo 7/ogo 237gg 217jogg 163gg 7jogg 217go" 

18 ^ 240 40 48 24 4 48 

_ ZM + _ 7jogo OlgogQ 7gggo y 

16 360 144 ^ 48 24 A’ 


(3.8) 


To avoid problems with the convergence of the series for the highest redshift 
objects, these relations are recast in terms of the new variable y = z/{l + z) 
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Hiy) = Ho 


2^1) y" 


1 + (1 + go) ?/ + (l + “ I' 

+ |(6 + 6go + 3 jo - 4gojo - 3g^ + 3g^ - sq) 2/^+ 

+ (l + 5'o ~ 2gojo + |<?o ~ I'^o) 1 /"^ + O [y^') 


dL { y ) = Ho ^ 


y + \{^-qo)y‘^ + \ (11 - 5go - jo + 3gg) y^+ 

(50 - 26go - 7jo + 21gg + lOgojo - 15go + sq) y"^ + O (y^) 


dA{y) = Ho^ 


1 / - 1 (1 + go) + 1 (1 - go + jo - 3gg) |/3+ 

+ ^ (“2 + 2go + jo — 3gQ + lOgojo — 15go + So) y"^ + O (y^) 


These formulae can become useful if one increases accuracy of the cosmolog¬ 
ical observations and includes datasets at higher redshifts. 

The formula (13.ip enables us to find the luminosity distance basing on 
given function H{z). Let us now solve the inverse problem and find the 
Hubble parameter as function of the luminosity distance [iQ] . In the spatially 
flat case differentiation of fl3.1l) results in 


djdi) ^ 1 + ^ dijz) ^ {7 + zf 

dz H{z) l + z djil + z)-dL 


where prime denotes derivative with respect to redshift Then 


q{z) 



-l + il + z)^ 


can be rewritten in the form 


(l + ^)^<(^) 

dj{l + z) -di 


q{N) 


H\N) 

H{N) 


d'i{N) + dj{N) 
dj{N) + dL{Ny 


Here primes denote derivatives with respect to iV = Ina = — ln(l -t- z). 
For arbitrary geometry 


1 d- UKdL{z)dj{z)/{l -F z) _ (1 -F zfd'j{z) 

l + ujKd\{z)/{l +zy d'^{z){l + z) - di 
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3.2 Source counts 


Let us cite without derivation a useful formula for the number of astronomical 
sources with redshifts in the range {z, z + dz) [201 02] 

-^0 

Here n is the number of sources in a unit proper volume. This equation is 
applicable to all Friedmann models. 


3.3 Horizons 


In the present section we analyze connections between the fundamental cos¬ 
mological parameters with dimension of length: the particle horizon Lp, the 
event horizon and the Hubble’s radius Rh- The reader is assumed to be 
acquainted with dehnition of these quantities. As we shall see, connections 
between these quantities depend on (and in considerable measure are defined 
by) the DP, i.e. the character of the cosmological expansion. Moreover, 
global picture of the observable Universe and its time evolution depend on 
the DP as well. Indeed, due to the Hubble’s law, the galaxies situated on 
the Hubble’s sphere recede with light speed. Velocity of the Hubble’s sphere 
equals to time derivative of the Hubble’s radius Rh = c/H, 



H 


d 

c— 

dt 



c /a \ 


c(l + ()). 


(3.9) 


As one can see from fl3.9ll the Hubble sphere contracts when g < — 1, remains 
stationary when g = — 1 and expands when g > — 1. We should stress the 
Hubble sphere in general does not coincide with the horizons, except when 
it becomes degenerate with the particle horizon at g = 1 and with event 
horizon at g = —1. 

Let us now analyze [19] kinematics of the Hubble’s sphere (which defines 
boundaries of the observable Universe) for different regimes of expansion of 
the Universe. 

In the Universe with decelerated expansion (g > 0) the Hubble’s sphere 
has velocity exceeding the light speed by the quantity eg and thus it overtakes 
the galaxies situated on its surface. Therefore the galaxies initially situated 
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outside the Hubble’s sphere will initially enter inside. Galaxies at distance 
R > Rh are later R < Rfj, and their superluminal recession in the course of 
time becomes subluminal. The light emitted toward the observer by a galaxy 
outside the Hubble sphere recedes until it enters inside the sphere. Therefter 
it starts approaching us and becomes available to observations. Therefore, 
all decelerating Universes lack event horizons unless they terminate at some 
future time. 

In uniformly expanding Universes g = 0 the Hubble surface and the 
galaxies situated on it have equal velocities. Thus number of galaxies in the 
observable Universe remains constant. Then both particle horizon and event 
horizon are absent in such Universes. 

In case of the accelerating expansion g < 0 the Hubble sphere has velocity 
which is less than the light speed by the quantity g and thus it falls behind the 
galaxies, and therefore number of them decreases inside the Hubble sphere. 
All accelerating Universes have the property that galaxies at distance R < 
Rh are later R > Rh, and their subluminal recession in the course of time 
becomes superluminal. Light emitted outside the Hubble sphere recedes from 
the observer and can never approach the observer. There are events that can 
never be observed, and such Universe have event horizons. 

The above mentioned situation generates the following question of prin¬ 
ciple. Should we treat as real the galaxies which were initially inside the 
Hubble sphere and then became unavailable to observation? It is the possi¬ 
bility to check a physical theory experimentally which distinguish the physics 
from metaphysics. Progress in the experimental technique pushes the bound¬ 
aries of physics, winning over more and more objects from the metaphysics. 
Atoms, elementary particles and black holes transformed from metaphysi¬ 
cal objects into physical ones during the last century. Being inhabitants of 
evidently expanding Universe, we face the opposite situation. As well as in 
decelerated expanding Universe, there are galaxies so distant that no signal 
from them can be presently detected by terrestrial observer. However if the 
cosmological expansion accelerates then we recede the galaxies with super¬ 
luminal velocity. Therefore if their light have not reach us till now then it 
will never do. Are we treat them as a physical object (because the galaxies 
are products of the Big Bang as well as ours) or a metaphysical one (due to 
impossibility to observe them)? Those who believe science hction to be the 
realization of unlimited fantasy are quite mistaken. Science hction is dull 
and lacks any hight of fantasy compared with cosmology. 

Let us now consider several examples which demonstrate kinematics of 
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the Hubble sphere in popular cosmological models. In the power-law models 
(a oc t”) we have 


H 


n 

i’ 


q = 


1 — n 


R 



Rh — 


c 

n 


In the matter-dominated Universe n = 2/3, q = 1/2 the Hubble sphere 
expands at velocity 3c/2 and overtakes the comoving galaxies at relative 
velocity c/2. In the radiation-dominated Universe n = l/2, g = l the Hubble 
sphere expands at velocity 2c and overtakes the comoving galaxies at relative 
velocity c. 

When the scale factor a{t) oc exp(hft), H = const (de Sitter Universe), 
g = — 1, the Hubble sphere has a constant radius. Galaxies cross the Hubble 
surface at velocity c. light emitted by these galaxies reaches the observer in 
the inhnite future with inhnite redshift. All events outside the Hubble sphere 
can never be observed, and the Hnbble sphere acts as an event horizon. 

Using the dehnitions of the particle horizon and the event horizon one 
hnds 


dLp d 
dt dt 

dLe d 
dt dt 



1 - 

a{t) 

J a{t') 

2 

0 


OO 

f dt' 

a{t) 

J ait') 

t 


— Lp{z)H{z) -|- 1; 
= L,{z)H{z) - 1. 


(3.10) 


One can see that observable part of the Universe expands faster than the 
Universe itself. In other words, the observed fraction of the Universe always 
increases. Indeed, the particle horizon at distance Lp recedes with velocity 
Lp, while the galaxies at the particle horizon recede at velocity HLp. Hence 
the horizon overtakes the galaxies with the speed of light c. For Universes of 
constant q 

Rh = 1 + 0') Lp = 1 -\—, —H = —. 

q Kh q 

When g < 0 no particle horizon exists. The light cone in this case extends to 
t = —oo. When 0 < g < 1 the Hnbble sphere lies inside the observe Universe 
Rh < Lp and bodies receding at velocity c at Hnbble sphere have hnite 
redshift. In the radiation-dominated Universe, g = 1, and the Hnbble sphere 
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and observable Universe have the same size. Differentiating the relations 
fid.iup w.r.t. time, one obtains 



(3.11) 


Let ns cite as an example the expressions for the particle horizon in the 
single-component Universe filled with non-relativistic matter [201 02] . In this 
case 


/c = 0, qo = 1/2; 
k = l, go > 1/2; 




arcsinh 

in 


0 


0 


k = -1, go < 1/2. 


\/i-2go 


V 


To conclnde we would like to stress one more time that interconnection be¬ 
tween the Hubble radius and the particle horizon determined by the expan¬ 
sion type (i.e. the DP) plays a principal role in understanding of casual 
connections between different regions of the Universe. We cite below an 
extremely bright discussion of this question in [19]: 

”How, in the Universe of age t can causally connected distances of L 3> ct? 
Let two comoving bodies be separated by a distance L sufficiently small that 
each lies in the observable Universe of the other. Each body remains there¬ 
after permanently in the other’s observable Universe, and the ratio L/{ct) 
during expansion depends on the behavior of the Hubble sphere. 

In a decelerating Universe the Hubble sphere expands faster than Uni¬ 
verse, and a body at distance L either is inside or will soon be inside the 
Hubble sphere. Hence, any two bodies must eventually recede from each 
other at subluminal velocity, and the ratio L/{ct) will then decrease in time. 
In an accelerating Universe the Hubble sphere expands slower than Universe, 
and a body at distance L either is outside or will soon be outside the Hubble 
sphere. Hence, any two bodies must eventually recede from each other at 
superluminal velocity and the ratio L/{ct) will then increase in time. 

How can causally connected distances of L S> ct exist? The answer is that 
the Universe passes through a period of accelerated expansion, and causal 
connections of L < ct, established before acceleration, expand superluminally 
outside the Hubble sphere... A period of accelerated expansion distend all 
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previously established causal connections and increases the distance to the 
particle horizon.” 
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Chapter 4 

The Effects of a Local 
Expansion of the Universe 


4.1 General Description 


One would expect cosmological expansion to have a significant effect on the 
dynamics of massive objects |l3l HU 051 06], 07]. Considering the radial 
motion of a test particle in a spatially-flat expanding Universe it is easy 
to show that in the Newtonian limit the radial force F per unit mass at a 
distance R from a point mass m is given by 

7T) 

F = (4.1) 

Thus, the force consists of the usual 1/R^ inwards component due to the 
central (point) mass m and a cosmological component proportional to R 
that is directed outwards (inwards) when the expansion of the universe is 
accelerating (decelerating).The latter formula has evident origin. In order to 
describe the cosmological expansion one commonly uses two sets of coordi¬ 
nates: the ’’physical” (or Euler) coordinates {R, 6, ip) and comoving (or hxed, 
Lagrangian) coordinates (r, 9, </9)0. The two sets are related by the formula 
R{t) = a{t)r. Therefore a point which is hxed w.r.t. cosmological expansion, 
i.e. with constant coordinates {r,d,ip), has additional radial acceleration 


fR 

dR 


expansion 



a 


-qH^R. 


(4.2) 


'^The angular coordinates are the same for both sets as the cosmological expansion is 
assumed to be radial 
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Let us dwell a bit on the simplest consequences of the additional acceleration 
appeared in fld.ip due to the expansion of Universe. Consider a Universe 
which contains no matter (or radiation), but only dark energy in the form of 
a non-zero cosmological constant A. In this case, the Hubble parameter and, 
hence, the DP become time-independent and are given hy H = a/A/ 3 and 
q = —1. Thus, the force fl4.ip also becomes time-independent, 

777 1 

f = + (4.3) 

For the case of spatially hnite (i.e. non-pointlike) spherically-symmetric mas¬ 
sive objects fl4.3p is replaced by 


F 


M{R) 




(4.4) 


where M{R) is the total mass of the object contained within the radius R. 
If the object has the radial density p{R) then 


R 

M{R) = j 4:7ir‘^p{r)dr. 

0 

Although the de Sitter background is not an accurate representation of our 
Universe, the SCM is dominated by dark-energy in a form consistent with a 
simple cosmological constant. Even in the simple Newtonian case fl4.4p . we 
see immediately that there is an obvious, but profound, difference between 
the cases A = 0 and A 7 ^ 0. In the former, the force on a constituent 
particle of a galaxy or cluster (say) is attractive for all values of R and tends 
gradually to zero as i? —>■ cx) (for any sensible radial density prohle). In the 
latter case, however, the force on a constituent particle (or equivalently its 
radial acceleration) vanishes at the hnite radius Rp which satishes 

Rf = [3M{Rf)/A]^^\ 


beyond which the net force becomes repulsive. This suggests that a non-zero 
A should set a maximum size, dependent on mass, for galaxies and clusters. 

From fl4.4p . in the Newtonian limit, the speed of a particle in a circular 
orbit of radius r is given by 

V(R) = (4,5) 
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from which it is clear that no circular orbit can exist beyond the radius Rp- 
Let us now consider |1B] a more realistic, time-dependent model for the 
background cosmological expansion, treating the central massive object as 
a point mass. In m there was obtained the metric for a point mass m 
embedded in an expanding cosmological background. In the spatially-flat 
case, which is a reasonable description of our Universe, and using physical 
coordinates, the metric takes on the form 


ds^ = 


2m r, r,, ^ 

1 - — - R^H\t) 


( 2m 

dt^R2RH{t) ( 1 - 


2 m \ 

~r) 


-1 


- Ri^dRr. 


(4.6) 


This metric is only applicable outside r > 2m. Nonetheless it is appropriate 
to use it for our region of interest r m. In 071 there was also obtained 
the force per unit mass required to keep a test particle at rest relative to the 
central mass m. 


^-RH\t) RH\t){q{t) + l)M 

F* = _ ^ _ — _^_V (4 7) 

(1 - f - 1/2 (1 - I - R^H'^{t)f/'^ 

In the region m -C r l/i7(t) taking into account that F = —F*, one hnds 
the radial force on a unit-mass test particle as simply that given in fld.ip . 

The cosmological force component —q{t)H‘^{t)R in the general case is 
time-dependent. For example, for the standard LCDM cosmology, the cos¬ 
mological force term reverses direction at about z = Zt ^ 0.7, changing from 
an inwards directed force at high redshift (decelerating expansion ) to an out¬ 
wards directed force at low redshift (accelerating expansion). Moreover, the 
dominance of the dark energy component will increase and so the expansion 
will tend to the de Sitter background model considered earlier, for which the 
cosmological force term is time-independent. 

The time-dependence of the cosmological force term in the general case 
leads to the result that the important structure parameter Rp is also time- 
dependent. For a central point mass, this is given explicitly by 


rp(t) ^ 


m 


1/3 


qm^it)\ 


(4.8) 
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provided the universal expansion is accelerating, so that q(t) is negative. 
Time-dependent cosmological force term will act essentially differently (de¬ 
pending on its sign) on the formation and structure of massive objects (galax¬ 
ies or galaxy clusters) as compared with the simple special case of a time- 
independent de Sitter background; in the latter case, the expression fl4.8p 
reduces simply to Rp = (3m/A)^/^. At low redshifts, where the dark energy 
component is dominant, we might expect that the values of Rp obtained 
using fl4.8p will not differ signihcantly from those obtained assuming a de 
Sitter background. Clearly, if the expansion is decelerating (forever) then 
the force due to the central mass m and the cosmological force are both 
directed inwards and so there is no radius at which the total force vanishes. 

The radius R = Rp does not necessarily corresponds to the maximum 
possible size of the galaxy or cluster: many of the gravitationally-bound par¬ 
ticles inside Rp may be in unstable circular orbits. Therefore it is important 
to know the so-called ’’outer” radius, which one may interpret as the max¬ 
imum size of the object, it is the one corresponding to the largest stable 
circular orbit Rs- The radius Rs may be determined as the minimum of 
the (time-dependent) effective potential for a test particle in orbit about the 
central mass [48] . 

Remaining within the frames of the Newtonian approximation (a weak 
gravitational held and low velocities), the equation of motion for the test 
particle reads 

(4.9) 


This is simply the Newtonian radial force expression fl4.ip with the inclusion 
of a centrifugal term. The radius of largest stable circular orbit in this case 


IS 


Rs{t) = 


m 


4q{t)H^t) 


1/3 


(4.10) 


Comparing this expression with that for Rp{t) in fl4.8p . we see that although 
both radii are time-dependent, they are related by a constant factor: at 
any epoch, the radius Rs(t) of the largest stable circular orbit lies a factor 
4^/^ 1.6 inside the radius Rp{t) at which the total radial force on a test 

particle is zero and the circular velocity vanishes. Using Rs{t) instead of 
Rp(t) enables us to correct our estimates of the maximum possible sizes of 
galaxies or galaxy clusters [46] . 

If we consider Rs{t) as the maximum possible size of a massive object 
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at cosmic time t, and assume that the object is spherically-symmetric and 
have constant density, then it follows from fId.lOp that there exists a time- 
dependent minimum density (due to maximum size) for objects, given by 


. , 3m 3q(t)H^(t) 

= 4 ^ = 


(4.11) 


As well as fld.lOp . the latter relation is valid only for q{t) < 0 (accelerating 
expansion). For q(t) > 0, Pminit) = 0. It is easy to see that 


Pmin (^) 
peril (^) 


= -8g(t). 


(4.12) 


For the current moment of time pminiio) ~ 4.4pcrit(^o)- One more im¬ 
portant characteristic - the minimum fractional density 5min = [Pminif) — 
Pmii)]/Pmit^ contrast - is closely related to the DP. As Pmii) = ^mif)peritif)^ 
then 


^min (^) 


1 + 


8q{t) 


(4.13) 


For the current moment of time (D^o ~ 0.3, go ~ —0.55) one hnds |16] that 
^minO ~ 14 . 


4.2 ”In an expanding Universe, what doesn’t 
expand?” 

The results obtained in the previous section formally related to the objects 
with cosmological spatial dimensions (such as galaxies and galaxy clusters). 
Because of importance of those results we would like to follow [HI ED] in 
detailed consideration of universal physical aspects of the influence the ex¬ 
pansion makes on the local physical objects, taking into account both the 
non-relativistic and relativistic effects. 

It is principal to answer the following question; if space itself is stretching, 
does this mean that everything in it is stretching? The traditional answer is 
that the ’’bounded” systems do not take part in the expansion. But if the 
whole space is stretching, how the bounded systems can avoid even minimal 
stretching? And what does it mean ”to be slightly affected” for a bounded 
system? Are the bounded system to be stretched with lower rate? Are the 
weakly bounded systems to be stretched stronger? 
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Price and Romano [H] made an attempt to answer those qnestions with 
the help of a simple model — a classical ’’atom”, composed of a negatively 
charged ’’electron” of negligible mass, orbiting aronnd a positively charged 
heavy ’’nnclens”. Let ns place the classical atom in a homogeneons Universe, 
where expansion is described by the scale factor a{t). Shall this expansion 
force the atom to grow in size, i.e. to increase the radins of the hrst Bohr 
orbit? Does the answer depend on size of the atom? 

For the concrete model nnder consideration the qnestion whether the 
atom joins the cosmological expansion can be reformnlated in the follow¬ 
ing way: what coordinate is hxed for the electron - the physical R or the 
comoving r? Or an intermediate case is realized? If R = const, then the 
cosmological expansion is absent for onr model atom. If r = const then the 
atom completely joins the cosmological expansion. 

Price have shown that the answer on the latter qnestion contains both 
expected and nnexpected featnres. The expected one: it is the relation 
between the kinematic characteristics of the expansion and the electrostatic 
interaction that determines which case (the hrst or the second) is realized. 
It tnrned ont that snfhciently weakly bonnded electron will expand together 
with the whole Universe, i.e. r = const, while more tightly bonnded one 
after initial perturbation of the orbit ignores further expansion and conserves 
R = const. It is unexpected that the intermediate case is absent. As Price 
notes, the ’’all-or-nothing” situation takes place. 

The model by Price places the massive nucleus at rest in the origin of 
the spherical coordinate frame {R,d,ip). Position of the electron of mass 
m orbiting in the equatorial plane 6 = 7r/2 is described by the functions 
R{t),ip{t). As there are only radial forces acting on the electron, then its 
angular momentum L = mR^ip is conserved, and we dehne the integral of 
motion for the unit mass electron 

L = R^p. (4.14) 

In the Newtonnian approximation in presence of the cosmological expansion 
a radial force per unit mass (teh equation of motion for R{t)) takes on the 
form 

= ( 4 - 15 ) 

Here C oc qQ (the proportionaly coefficient depends on the choice of the unit 
system, for example in SI units C = qQ/{A-KSom) . Of course the equation 
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fl4.15p up to redefinition of the constant C coincides with the equation fl4.9p 
for a test particle in the gravity held of a galaxy. 

The comparative strengths of the electrostatic and cosmological terms in 
Eq. f|4.15p can be usefully cast as a comparison of time scales for atomic and 
expansion effects. Let us dehne a characteristic atomic time scale tatom as a 
combination of the parameters (L and C) relevant to the electron’s motion, 
tatom = L^/C^- In the absence of expansion effects the time for the electron 
to complete a circular orbit is 27itatom ■ 

Note that the dynamics of the expansion is coded in functional form of 
a{t). Choice of this function means hoice of the expansion kinematics, dic¬ 
tated by the dynamics. The concrete choice of a{t) is irrelevant to answer our 
question of interest, so choosing a{t) one can use the convenience arguments 
if remaining in frames of realistic expansion laws. So let us consider the sim¬ 
plest type of the accelerated cosmological expansion — de Sitter cosmology. 
In this case a{t) = exp(t/f//), q = —1, H = const = t~^ and the equation 
fl4.15p takes on the form 


i?2 + i2^ + ^3- 


The hrst integral of the latter reads 


r . 1 ,^2 


C 

R 




(4.16) 


(4.17) 


The electron then feels the effective one-dimensional potential 


V{R) 


L2 c r^ 

^~R~ 


(4.18) 


Graphical analysis of the potential fl4.18p at different values of the parameters 
tatom and tn allows to answer the above posed question about the charac¬ 
ter of the electron’s motion. A graph of the potential for various values 
T = tatom/tn is showu in Fig J4.ll where the dimensionless potential VLi^ 
versus the dimensionless radial distance CR/L? is shown. Each curve is la¬ 
beled with the value of the parameter r that determines how strongly the 
cosmological expansion affects the evolution of the atom. The larger the 
value of r, the larger is the effect of expansion. 

Let us analyze behavior of the electron for different values of the parame¬ 
ter T. Expansion is absent for the top curve, for which r = 0. In this case, the 
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Figure 4.1: Effective potential for exponential expansion. Curves are marked 
by the value of the parameter r = tatom/tn- The curve labeled 0 is the no¬ 
expansion potential, for which r = 0. The dashed line shows the alignment 
of the minimum for the no-expansion potential with the local maximum of 
the potential for r = 0.25 [44] . 

electron is always trapped in the potential well, i.e., it is permanently bound. 
If the electron begins at the bottom of the well {R = I? jC^ E = —C^/(2L^)) 
it will remain in a circular orbit at that radius for all time. For any larger 
value of E, the electron will orbit in an ellipse. For nonzero values of r, the 
potential at large R eventually becomes negative and decreasing, thus repre¬ 
senting a dominant outward force. Consequently, an electron at a sufficiently 
large distance from the nucleus will be driven to an even larger distance. The 
important question is whether the electron will ever get to this region of dom¬ 
inant outward force. The answer is contained in the shapes of the curves in 
Fig J4.ll We hrst consider the situation where the electron sits at the bottom 
of the no-expansion potential well and is ’’surprised” when the expansion 
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is suddenly turned on. Thus, the electron has energy E = —C‘^/{2L‘^) and 
hnds itself under the influence of one of the expansion potential curves with 
r > 0. In this scenario, there is a critical value of r = 0.25, above which 
the electron will be accelerated outward by the cosmological expansion. As 
shown by the dashed line in Fig i4.ll this critical value occurs when the local 
peak in an expansion potential has the same value as the lowest point in the 
no-expansion well. For 0 < r < 0.25, the electron will remain trapped in an 
approximately elliptical orbit. 

A different scenario can also be envisioned. Imagine the electron is sitting 
at the bottom of an expansion potential well. In this case, the electron will 
remain at a fixed R (the bottom of the well), assuming such a local minimum 
actually exists. However, as shown in Fig J4.ll there is a critical curve that 
separates potentials that have a local minimum from those that do not. This 
curve has r = 3\/3/16 ~ 0.324. This value can be easily obtained using the 
following arguments. 

Extrema of the effective potential in which the test particle moves occur 
at the i?-values for which (PR/dR? = 0, namely the solutions of 


C R 


0 . 


Consider the function 


y = -CR + L^. 


(4.19) 


(4.20) 


This polynomial (naturally of i? > 0) has a unique extremum — a minimum 
at 



The existence condition for real roots of the equation (I4.19p reads y{R*) < 0. 
The critical value of the parameter r at which the minimum of the effective 
potential fld.lSp disappears, corresponds to the condition y{R*) = 0 and then 
T* = 3^/3/16 0.324. 

Of course the redehnition m = C and H = map all the above obtained 
results on the case of the test particle in the gravity field of spherically 
symmetric galaxy considered in the previous section. 

Such a qualitative analysis allows us to understand why the atom has 
an all-or-nothing behavior. The electron either is, or is not, trapped in the 
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potential well. Correspondingly, the atom either expands or does not; there 
is no ’’partial expansion” possible. Underlying this graphical understanding 
is a broader but less precise heuristic explanation of the all-or-nothing effect, 
an explanation that applies regardless of the specihc form of the expansion. 
The cosmological expansion term 

-R{-qH^R) 

a 

increases at large physical distances R from the nucleus, whereas the cen¬ 
trifugal and electrical forces both decrease. This implies a sort of instability 
with respect to expansion. If the electron moves suf?ciently far from the 
nucleus, the expansion term becomes more important and this pushes the 
electron even further away. 

We can get yet another viewpoint on the bound/unbound issue by nu¬ 
merically solving Eq. fl4.16p . If we start the computation with the electron 
at the bottom of an expansion well, the results are in agreement with the 
predictions of the analysis based on Fig M.ll — the electron remains at hxed R. 
More interesting is the ’’surprised electron” scenario discussed above (with 
E = —C/{2L‘^)). The results, shown in Fig l4.21 are in accord with the anal¬ 
ysis based on Fig l4.ll For r slightly greater than the 0.25 critical value, the 
physical radius R of the atom grows exponentially after an initial hesitation. 
In contrast, for r slightly less than this critical value, the electron remains 
trapped in an approximately elliptical orbit and is not dramatically affected 
by the exponential expansion. 

4.3 Tethered Galaxy 

As we have seen above, the unstable character of the cosmological accel¬ 
eration leads to different observable effects. They can be divided into two 
groups. The hrst is due to variation of relative distance between the mem¬ 
bers of the cosmological expansion regardless the possibility to register this 
variation. The second is connected with analysis of the additional red shift 
due to the relative motion. The DP plays a crucial role in both cases. In¬ 
stead of scholastic disputes on the nature of the red shift which is connected 
with numerous misunderstandings |51], we consider a simple model called 
the ’’tethered galaxy” [52] . 

Let a distant galaxy in the expanding Universe recedes with velocity V 
given according to the Hubble law by U = HR. Imagine that we separated 
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Figure 4.2: Radial coordinates as a function of time for exponential expan¬ 
sion. On the left is the case for r = 0.2505 for which the electron’s comoving 
radius r remains approximately constant after an initial decrease to about 
2% of its initial value. Due to the exponential increase in a{t), the physi¬ 
cal radius R grows without bound. On the right is the radial kinetics for a 
slightly smaller value of r = 0.2495. In this case, the electron remains bound 
in an approximately elliptical orbit with the physical radius oscillating be¬ 
tween values near the original atomic radius. The coordinate radius R in this 
case falls off exponentially [H] . 
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a small test galaxy from the Hubble flow and tied it to the observer galaxy 
so that the physical (proper) distance between them remains constant. We 
will not touch the practical realization of this project, as it always can be 
thought as we apply to the tethered (or now not tied) galaxy a velocity 
directed toward the observer so that this velocity precisely compensates the 
bubble velocity of the expansion. After we remove the tie (or switch off the 
jet engine) there is set the initial condition R = 0. We dehne total velocity 
of the test galaxy as the time derivative of the proper distance 

R^ R hr “I" QjT “1“ f^ec* 


The peculiar velocity Vp^c is the relative velocity w.r.t. the comoving 
coordinate frame which the test galaxy is extracted from. This velocity 
corresponds to our common understanding of velocity and it must be less 
than the light speed. It is presence of the peculiar velocity which makes the 
derivative r 7 ^ 0. The recession velocity Kec of the test galaxy equals to 
the velocity of the Hubble flow on the distance R and it can be arbitrarily 
high. We will assume below that the test galaxy has negligibly small mass. 
According to the construction, the tethered galaxy initially has zero total 
velocity Rq, or 

Vrec ~ Vpec^ ^ hgTo = CIqTq. 

With these initial conditions we untie the galaxy and let it sail freely. 
What will it do then; approach to the observer (in the origin of the coordinate 
frame), recede or keep the distance? 

It is principal that the momentum p w.r.t. the comoving reference frame 
decays as 1/a. Let us consider the non-relativistic case for simplicity ( [S2] 
contains both non-relativistic and relativistic cases). For the non-relativism 
p = mVpec and, therefore. 


^pec h/eco/®5 



In order to answer the above posed question one should analyze the time 
dependence of the proper distance R between the test galaxy and the observer 


dt 


R = aro 1 — ho 


(4.21) 
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Figure 4.3: Graphical interpretation of the tethered galaxy problem. 


The integral fl4.2ip can be evaluated numerically using dt = da/d and ho 
obtaijned immediately from the hrst Friedmann equation with hxed energy 
content of the Universe. Figure 14.31 demonstrates solution for four differ¬ 
ent models: SCM = (0,3; 0,7), empty Universe = (0;0), 

Einstein-de Sitter model (r2m,GA) = (1,0) and (r2m,UA) = (0,3;0). Dif¬ 
ferent behavior in the considered models is determined by the composition 
of the Universe in each model. It is the composition of the Universe which 
determines the kinematics of the expansion (accelerated, decelerated or uni¬ 
form). For the four cosmological models we untether a galaxy at a distance 
of Rq = lOOMpc with total initial velocity equal to zero and plot its path. 
In each case the peculiar velocity decays as 1/a. Its final position depends 
on the model. In the SCM , accelerating Universe, the untethered galaxy re¬ 
cedes from us as it joins the Hubble flow, while in the decelerating examples 
(Dm, Ha) = (1;0) and (Dm,HA) = (0,3;0) the untethered galaxy approaches 
us, passes through our position, and joins the Hubble flow in the opposite 
side of the sky. In the (Hm, Ha) = (0; 0) model the galaxy experiences no 
acceleration and stays at a constant proper distance as it joins the Hubble 
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flow. 

We have considered above the behavior of the test galaxy in the Universe 
with hxed DP. The models where the DP can change sign during the time 
evolutions appear more realistic. Let us make some preliminary notions. 

It is not hard to see that the test (untethered) galaxy asymptotically 
joins the Hubble flow in any cosmological model which is always expanding. 
Joining of the Hubble flow takes place because 

Vtot = Vrec + = Kec + (4.22) 

a 

At a —)• cx) Vtot = R = Vrec = HR, and it is purely Hubble flow. Note 
that the galaxies join the bubble flow due to growth of the scale factor (or 
expansion of the Universe). Differentiating fl4.22p w.r.t.time one obtains 

R= {dr + dr) — j _ Vpec- = {dr + dr) — ar- = d. (4.23) 

a a a a 

Using the dehnition of the DP one hnds 

ii = -qH^R. (4.24) 

An alternative way to obtain fl4.24l) is to differentiate the Hubble law. This 
approach ignores the peculiar velocity and therefore it does not include ex¬ 
plicit cancelation of the two terms in (I4.24p . which represents more general 
consideration. The fact that both results are the same, shows that the ac¬ 
celeration of the test galaxy equals that of the comoving galaxy and there 
is now additional acceleration of the galaxy getting captured by the Hubble 
flow. 

Recall that i?o = 0 according to the chosen initial conditions. Therefore 
it is sign of acceleration R of the galaxy which determines whether the galaxy 
approaches us or recedes. The sign in its turn is determined by the sign of 
the DP according to fl4.24l) . Equation fl4.24p shows that the sign of R is 
determined by the sign of the DP q. If the Universe expansion is accelerated 
{q < 0) then R > 0 and the test galaxy recedes from us. If the Universe 
expansion is decelerated {q > 0) then i? > 0 and the test galaxy approaches 
us. At last in the case q = 0 the proper (physical) distance does not change 
as the test galaxy joins the Hubble flow. The consequence is that in the 
Universe where the expansion regimes interchange the time dependence of 
the distance between the test galaxy and the observer is very complicated 

(FigJiaD. 
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Figure 4.4: |52] Upper panels: = (0,95; 0,05). This particular 

model shows the effect of q (right panel) on the position of the test galaxy 
(left panel). Initially q > 0 and the proper distance to the untethered galaxy 
decreases. But q subsequently evolves and becomes negative, reflecting the 
fact that the cosmological constant begins to dominate the dynamics of the 
Universe. With g < 0 the acceleration R changes sign. This makes the 
approaching galaxy slow down, stop, and eventually recede. The dotted lines 
are fixed comoving coordinates. Lower panels: (f2m,f^A) = (2;0). Universe 
expands and then recollapses (a changes sign). 
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4.4 Inertia forces in the accelerated expand¬ 
ing Universe 

As we have seen above fl4.24p . as the Universe expands, the relative acceler¬ 
ation between two points separated by a distance R is given by —qH'^R. If 
there is a particle with mass m at each of the points, an observer at one of 
the masses will measnre an inertial force on the other mass of 

Finert = -mqH'^R. (4.25) 

Let ns estimate magnitnde of the force fl4.25p |53j . In the Universe filled 
by non-relativistic matter with zero pressnre and dark energy with the state 
eqnation p = WdeP the eqnation fl4.25p (in cnrrent moment of time) can be 
rewritten as 

Finert ^ mqHQR{l + 3wde^de)- (4.26) 

For WdC ~ -1, Qde ~ 0.7 

Finert = (3 X 10“^® sec“^)mi?. (4.27) 

For a first glance it seems that one can nse the obtained estimate for Finert 
to establish limitations for the state eqnation parameter for the dark energy 
[53] . The stability condition for a galaxy clnster in the expanding Universe 
can be rewritten in the form 


F 


grav 


> F, 


inert' 


(4.28) 


The inertial force fl4.25p can be rewritten as Finert = —{mRiF/6){p + 3p), 
= SttG. Accordingly, the Newton gravity which the point particle with 
mass m snffers is given by 


Fgrav G 


mM AttG 


TflRpidust- 


(4.29) 


It follows then that the stability condition for the galaxy cluster reads 


Fgrav > Finert- (4.30) 

Since pciust 200 p, we hnd, that stability of the cluster is provided by the 
condition Wde > — 70, which does not pose any noticeable constraint on the 
state equation parameter for the dark energy. 
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